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Abstract 

Using molecular dynamics computer simulations we study the dynamics of 
a molecular liquid by means of a general class of time-dependent correlators 
Sw(q,t) which explicitly involve translational (TDOF) and orientational de- 
grees of freedom (ODOF). The system is composed of rigid, linear molecules 
with Lennard-Jones interactions. The (/-dependence of the static correlators 
Sjp(q) strongly depend on I, I' and m. The time dependent correlators are cal- 
culated for 1 = 1'. A thorough test of the predictions of mode coupling theory 
(MCT) is performed for Sff(q, t) and its self part S^ )m (q, t), for I = 1, . . . , 6. 
We find a clear signature for the existence of a single temperature T c , at which 
the dynamics changes significantly. The first scaling law of MCT, which in- 
volves the critical correlator G(t), holds for I > 2, but no critical law is 
observed. Since this is true for the same exponent parameter A as obtained 
for the TDOF, we obtain a consistent description of both, the TDOF and 
ODOF, with the exception of I = 1. This different behavior for I ^ 1 and 
I = 1 can also be seen from the corresponding susceptibilities (x")?l(q,u)) 
which exhibit a minimum at about the same frequency u m in for all q and all 
/ ^ 1, in contrast to (x")u(q,io) for which oj' min 10 0J m in ■ The asymp- 
totic regime, for which the first scaling law holds, shrinks with increasing I. 
The second scaling law of MCT (time-temperature superposition principle) 
is reasonably fulfilled for I ^ 1 but not for 1 = 1. Furthermore we show that 
the q- and (/, m)-dependence of the self part approximately factorizes, i.e. 
Sli )m (q,t) Cl s) (t) F s (q,t) for all m. 

PACS numbers: 61.43.Fs, 61.20.Ja, 02.70.Ns, 64.70.Pf 
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I. INTRODUCTION 



In the last few years quite a few papers were published in which computer simulations 
were used to study the time dependence of the translational degrees of freedom (TDOF) in 
supercooled liquids. On the other hand, the orientational degrees of freedom (ODOF) were 
so far investigated in much less detail since the simulation and data analysis of systems in 
which the particles are molecules are quite a bit more involved than the ones in which the 
particles have no structure. However, since most real materials are of molecular nature and 
since experimental methods such as light scattering or dielectric measurements probe also 
the ODOF, it is important to understand how the dynamics of the TDOF and the ODOF 
are related to each other. Only by understanding this relationship it will be possible to 
make a correct interpretation of the experimental measurements and to gain insight into the 
nature of the glass transition, i.e. the dramatic slowing down of the relaxation dynamics of 
supercooled liquids upon approaching the glass transition temperature. A more thorough 
discussion of these connections can be found, e.g., in Ref. [[I]], were we also review some of 
the other work in this field. 

Very recently we have carried out a molecular dynamics computer simulation of a simple 
molecular system in order to make a detailed comparison between the dynamics of the 
TDOF and the ODOF [JTJ] . Each molecule in this system is dumb-bell shaped and consists 
of two Lennard- Jones particles that are separated by a fixed distance d. More details on the 
system and the simulation can be found in Ref. Q. In that paper we studied the time and 
temperature dependence of the orientational correlation functions 

^) = ^E(Wn(o)-M*))> , f>i , (i) 

n,n' 

and the self part C$ (t). Here u n (t) is the unit vector pointing along the molecular symmetry 
axis of molecule n and Pi(x) is the /-th Legendre polynomial. The relevance of these type 
of correlation functions is given by the fact that they can be measured in experiments. The 
main results of that paper were that the temperature dependence of the relaxation times 

(s) 

of Ci, C\ and the diffusion constant D were given by a power law with the same critical 
temperature T c but with critical exponents that depend on the observable. In addition we 
showed that the so-called time temperature superposition principle works well for C\ s \ if 
/ > 2. Thus we concluded that many of the predictions of mode-coupling theory (MCT) 
hold for these correlation functions, although certain discrepancies are present. 

In the preceding paper, subsequently called KKSI, we have investigated the time and 
temperature dependence of the translational degrees of freedom by studying quantities like 
the van Hove correlation function G(r, t) and the intermediate scattering function F(q, t) [|J. 
The main conclusion of that paper was that MCT is able to give also a good description for 
the time and temperature dependence of these correlation functions. 

As we will demonstrate below, the intermediate scattering function F(q, t) and the orien- 
tational correlation functions Ci(t) are just a special case of a more general type of correlation 
function, which involves the translational as well as the orientational degrees of freedom at 
finite wave-vector q, i.e. |g| > 0. The goal of the present paper is therefore to investigate 
the time and temperature dependence of these more general correlation functions, since it 
is these correlators which are needed for a more detailed description of the dynamics of a 
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molecular system. In addition these correlation functions can also be calculated directly 
within the framework of MCT (although such a calculation might be in practice quite in- 
volved) thus allowing to perform a more stringent test of whether MCT is able to give a 
correct description of the dynamics of the system investigated. 

Our paper is organized as follows: In the next section we will introduce the mentioned 
generalized correlation functions and will discuss some of their properties. Section [III 
presents the results and the MCT-analysis and the final section contains a summary and 
our main conclusions. 



II. CORRELATION FUNCTIONS 

We introduce a set of correlators which involves the one-particle density (including the 
angular dependence) for a molecular liquid of rigid, axially symmetric molecules: 

N 

P (x, n,f) = £ - 5{n, n n (t)) (2) 

71=1 

where x n {t) and O n (t) = {9 n (t), <f> n {t)) denote the center of mass position and the orientation 
of the n-th molecule at time t, respectively. Due to the non-Euclidean metric for the angles 
6 and 0, one must use the invariant delta function 5(Q,Q'). For this and other details of 
the theoretical description of molecular liquids, the reader is referred to the textbook by 
Gray and Gubbins ||. Expansion of p(x, fl, t) with respect to a product of plane waves and 
spherical harmonics Yi m (Q) leads to the tensorial density modes 

N 

Plm (q, t) = l 1 ^ Ylm{ttn{t)) , (3) 

n=l 

where I = 0, 1, 2, . . . and —l<m<l. The factor v4tt is used so that poo(q, t) equals the 
definition of p(q,t) for simple liquids and i l is introduced for convenience (see below). The 
corresponding correlators 

Si m ,i'm'{q,t) = j^{8p*i m {q,t) 8pi, m ,(q,0)) (4) 

of the fluctuation Spi m (q,t) = pi m (q,t) — (pi m {q,t)) vanish for (q,l,m) = (0,0,0), and are 
otherwise given by: 

A/jr 

Si m ,i>m>(q,t) = — i 1 '- 1 £ (exp \-iq- (x n (t) - x n ,(0))\ Y*jn n (t)) ^ m ,(fi n ,(0))) (5) 

n,n' 

which shows the explicit dependence on both, the TDOF and the ODOF. Its corresponding 
self part S[^ Vml (q, t) is obvious. 

Taking into account that loo = l/V^vr one obtains from Eq. @: 

Soo,oo(<f, t) 
^oo,oo{q) 
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i.e. the normalized density correlator for the center of mass positions, which was studied in 
KKSI. On the other hand, we find from Eq. ([|) for q = 0: 

Si mi i' m '(0,t) = Ci(t)S mm /Sw . (7) 

Here the addition theorem for the spherical harmonics || and the isotropy have been used. 
As already mentioned in the Introduction, this special case was investigated in Ref. [|TJ. 
Eqs. (H) and (|7|) hold for the corresponding self part, as well. 

Although it is not obvious how these correlators for 1,1' 7^ can be measured in real 
experiments for q 7^ 0, they are the basic quantities which enter the MCT for a molecule 
in a simple liquid || and for molecular liquids [f^-Hfl. To our knowledge, there exists only 



one computer simulation which considers g-dependent orientational correlators [12]. But the 



experimental relevance of these correlators considered in Ref. |L2| is unclear. The correlators 
given in Eq. (^) simplify a bit, if one uses the g-frame ||, i.e. q — % = (0, 0, q). In that case 
one obtains ||: 

Sim,i> m i(qo,t) = S^,{q,t) 8 mm i , (8) 

which differ from zero only for < |m| < min(/,/') . Since S{p(q,t) = S^, m (q,t), one can 
restrict oneself to m > 0. The introduction of i l in Eq. (H) makes Slft(q,t) a real quantity. 
The same properties hold for the self part as well. In the following we will present all results 
in the g-frame. 

Some of the equations that we will subsequently make use of have been given in KKSI 
and are not reproduced here. We will refer to the nth equation of that paper by (1-n). 



III. RESULTS 

This section is subdivided into two parts. The first part contains the results for the 
static correlators Slp(q), and the second one presents the dynamical correlators Slp(q, t) and 
Sft,(q,t). In the following we restrict the values of I and /' to 0, 1 and 2. 



A. Static properties 

The static correlators are shown in Figs. [I] - 13| for the lowest investigated temperature 
T = 0.477. First of all, it becomes obvious from these figures that S$(q = 0) is m- 
independent and diagonal in / and as it should be due to isotropy. A comparison of the 
various diagonal correlators in Figs. [I] and ^] with each other shows, that the correlators 
Si t (q) for I = 1 and 2 possess a significant g-dependence similar to that of S(q) = F(q, 0), 
in contrast to those for m^0. The same behavior was found for a system of dipolar hard 
spheres ||, although for that system the most prominent peak occurs for Sji(q) at q = 0. 
In contrast to 5(g) and ^(g), the correlator S^iq) has a rather broad maximum at g = 
with a height which is comparable to that at q' max ~ 7.3, the location of the main peak in 
5*22 (g)- In Fig. |3| we present the non-diagonal correlators S™,{q) with I ^ I'. First of all one 
recognizes that Sq 2 (q) is much larger than S^q) and S'^(q). This can easily be understood. 
If the molecules had "head-tail" -symmetry, then it can be shown that S]Ji(q) = 0, for 1,1' 
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such that I + V is odd. Since for our molecules this symmetry is only slightly broken, we 
expect S™,{c[) to be much smaller for I + V odd than for I + V even. 

The second point one recognizes from this figure is that the non-diagonal correlators 
SJp(q) can have the same magnitude than the diagonal ones. Hence, there is no reason why 
the former should be neglected in analytical calculations. For example, since the solutions 
of the MCT-equations for the time-dependent correlators S$(q,t) are determined by the 
static correlators S^(q), it might not be a good approximation to consider I = I', only. 



B. Dynamical properties 

We have investigated both, the self correlators for I — I' — 0,1, ... ,6 and the collective 
correlators for I = I' = 0, 1 and 2. Let us start with the self part S^ m (g, t). Often it is 



assumed (see e.g. Ul3| ) that the q- and (I, m)-dependence (where I = V) factorizes, i.e.: 

Sl° )m (q,t)=Cl s \t) F s (q,t) (9) 

with cf°(t) the self part of Eq. <g) and F s (q,t) = S { l )0 (q,t), the self part of Eq. (§). The 
reader should note, that Eq. (||) is assumed to hold for all m, and that the factorization 
is trivial for q = 0. To check the validity of Eq. (0) for q > 0, we show Sl^ m (q,t) and 
C z (s) (t) • F s {q,t) in Fig. | (7 = 1) and Fig. § (I = 2) for three different q- values and T = 
0.477. Although the factorization becomes worse with increasing q, it is still a reasonable 
approximation, even for q = 10.6. Furthermore, the quality of the factorization is better in 
the /3-relaxation than in the a-relaxation regime (at least for 1 = 2), and it also becomes 
better with increasing temperature. 

This approximate factorization does not necessarily mean that the coupling between the 
TDOF and ODOF is very weak. The comparison of Cl s \t) with F s (q,t) in Fig. | and 
Fig. |5| reveals the reason why S^ m (q,t) can be approximately factorized. For instance, 
C[ s \t) has decayed to 0.1 for t = 2 • 10 4 , whereas at this time the value of F s (q = 2.8, t) 
is still around 0.85, i.e. the ODOF relax much faster than the TDOF. This is consistent 
with our observation that in the time span of the orientational correlation time, as deduced 
from c[ s \t) at the lowest temperature, the average center of mass positions change only a 
fraction (about 30 %) of the mean distance between the molecular centers. We stress that 
this is different to the MD-simulation of supercooled water. There, F s (q,t) and Ci(t) relax 
on approximately the same time scale jn] . 



We now turn to the test of the various MCT-predictions (see KKSI). We find |1| that 
S[i(q,t) do not obey the second scaling law, i.e. the time-temperature superposition 
principle. This observation has already been made for the case q = 0, which shows that 
this type of correlation function does not follow the predictions of MCT. This situation 
is different for the correlation function Si^ m (q,t) with / > 2 for which the second scaling 
law holds reasonably well. The critical exponents (which are practically q- independent) for 
the divergence of the relaxation time, '-/{^ and 72^, is 1.8 and 2.45, respectively, where the 
latter value is fairly close to the one found for the TDOF, F s (q,t), which was 2.56 ||. The 
exceptional role for the correlators with I = 1 is due to the existence of 180°-jumps of the 
molecular axis , since the Legendre polynomial Pi (cos 6) is sensitive on reorientations by 
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180°. The same is true for all P/(cos#) with I odd. But the weight of Pi(cos8) for 9 « 0° and 
6 ~ 180° decreases with increasing I. Since the second scaling law holds for I ^ 1, we can 
restrict ourselves in the following to the analysis of the correlation functions at the lowest 
temperature. 

In Fig. ^| we investigate the validity of the first scaling law [Eq. (1-4)]. This is done for 
q = by fitting Ci S \t) with the critical correlator G(t). We remind the reader that this 
fit is performed for fixed values A = 0.76 and t a = 69 as obtained from the similar fit of 
F(q max ,t). More details on this analysis can be found in section [TV] of the preceding paper 
J|]. For I > 2 (Fig. §) the critical correlator fits the data very well over about two decades 
in time. This range, however, becomes smaller with increasing I, which may indicate that 
corrections to the asymptotic law become more important for large /. If one uses A and t a 
(cf. (1-4)) as free fit parameters, the resulting fits follow the data longer by additional one 
to two orders of magnitude in time. (We note that even C[ (t) can be fitted reasonably 
well with G(t). Since we have shown in Ref. [|lj that for this correlation function the first 

(s) 

scaling law does not hold, one might argue that it does not make sense to analyze C\ in 
the way proposed by MCT. However, we find that the violation of the second scaling law is 
only weak and therefore it is not unreasonable to make such an analysis.) The so obtained 
values for A increase towards one with increasing I and reach, e.g., 0.97 for I = 6. We also 
mention that we do not observe a critical law, Eq. (1-6), the reason for which is likely the 
strong influence of the microscopic dynamics on the early /^-relaxation regime. 

We have found that these results do not change significantly for Sip m (q, t) if q > 0. From 
the fit with von Schweidler law plus corrections, Eq. (1-9), (not shown in Fig. |]) one can 
deduce the critical nonergodicity parameter //* (?), the critical amplitude h^ m (q) and 
the correction h\i m (q) which are shown in Fig. [7| for I = 1, 2 and 6, for the case m = 
(see KKSI for the difference between (h(q), h^ 2 \q)) and (h(q), h^(q))). We note that the 
result for I = 1 was obtained for A = 0.76 and a shift of the time scale to t' a = 10. Due 
to the approximate factorization property, the g-dependence of //j (g) is given by that of 
f( s ' c \q) = foo (q). The functions fu' c ^ m (q) decrease with increasing /, as expected from 
Fig. || The variation of the critical amplitude h^ m (q) and the correction m (q) with q is 
similar to that for I = V = (cf. Fig. 13 of KKSI) with the exception that these quantities 
do not vanish for q — > 0. 

The a-, (3- and the microscopic time scale can be better visualized from the imaginary 
part (x^ )/?(?) LJ ) °f the dynamical susceptibility as a function of uj, which is shown for 
m = in Fig. ^ for q = q max , I = and q = 0, I = 1, 2. The microscopic peak is at about 
u — 1 for all these values of I. Whereas the position of the a-peak and the location of the 
minimum (for low temperatures) are approximately the same for I = and 1 = 2, these 
positions are shifted to higher frequencies by about one decade for 1 = 1. We believe that 
this shift relates to the 180°-jumps of the molecules (see Ref. because these jumps do 
not affect the correlators with even I, but those with odd value of /, and particularly those 
with 1 = 1. 

The rest of this section is devoted to the discussion of the collective correlators Sjf (q, t), 
which are presented in Fig. |9] for q = 2.8 (the position of the main peak of S^q) (cf. Fig. [I])) 
and in Fig. |TTJ for g = 6.5 (the location of the main peak of S(q) = S® (q) (cf. Fig. [I])). Note, 
that, due to symmetry (cf. section ||), there are only two and three independent correlators 
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for I = 1 and / = 2, respectively. These correlators exhibit a strong m- dependence, in 
contrast to slp m (q,t). The reader should also note that ^(g, i) < ^(g, £) for q = 2.8, 
whereas S^q, t) > S^q, t) for q = 6.5. These inequalities are related to the fact that 
S^q) has its main peak at q = 2.8 where S^g) does not have a maximum, whereas (q) 
has its main peak at q = 6.5, where S^q) is close to a minimum. Similar considerations 
hold for the m- and g-dependence of S^(q,t). These observations make it obvious that a 
factorization [cf. Eq. (|9|)] does not work for the collective correlators. 

The test of the second scaling law is shown in Fig. [Ll| for q = 2.8, m = and I = 1,2. 
As already found for C[ • (t) and C/(i), i.e. the correlation functions for q = 0, this scaling 
law holds for I = 2 but not for 1 = 1. We define the a-relaxation time Ti mtq (T) as the time 
it takes S™ (g, r/ m>? ) to decay to the value of 1/e. The temperature dependence of Ti m>q {T) 
is shown in Fig. |T3. Fixing T c = 0.475, the a-relaxation times obey a power law (1-10) over 
about 2-3 decades in time. For the corresponding exponent 7 one obtains approximately 
1.9 for I = 1 and 2.5 for / = 2 with no significant g-dependence. Again the 7-values for I = 2 
(and the same remains true for I = 3, . . . , 6) fit with that for I = 0, which was around 2.55 
(see KKSI), whereas the value of 7 for / = 1 is quite different. 

The test of the first scaling law by fitting the time dependence of S]£(q,t) with the 
critical correlator is done in Fig. [TS| for I — 2, m — 0. This fit (again with A = 0.76 
and t a = 69) works well for different values of q. From the fit with the von Schweidler 
law plus correction, Eq. (1-9), (not shown in Fig. [13]) we compute the critical nonergodicity 
parameter fu' m (q), the critical amplitude h%(q) and the correction h^ m (q) } shown in Figs. |14 
and [15] for, respectively, / = 1 and 1 = 2. Although we have seen, that / = 1 is rather 
special, we have analyzed the corresponding correlators at the lowest temperature and have 
included its result. For reference we also show in Figs. [13] and |15] the static correlator 
S'a(q) and the a-relaxation time Ti mtg (T) for T = 0.477. These quantities possess the same 
characteristic g-dependence already found for the corresponding quantities of the TDOF, 
i.e. for l = l' = m = m' = (cf. Figs. 18 and 19 of KKSI). This means that (i) T\ mA 
and fii m {q) are in phase and h^{q) and h^ m (q) are in anti-phase with S™ (g) and (ii) the 
correction hff m (q) is smallest at that g where S™(q) has its main peak. This latter fact is 
well pronounced for (l,m) = (1,0) and (l,m) = (2,0) and less for the others, because there 
also the g-dependence of S]^(q) is less pronounced. 



IV. DISCUSSION AND CONCLUSIONS 

For a system of diatomic and rigid molecules interacting via Lennard- Jones potentials we 
have investigated by means of a MD-simulation the time and temperature dependence of a 
general class of g-, (l,m)- and (/', m')-dependent correlators. These correlators Sim^m/^t) 
contain the TDOF and ODOF explicitly. 

The static correlators S]]}(q) in the g-frame are not diagonal in / and /'. Whereas those 
with / + /' odd are smaller than S(q) = SQ (q) by about one order of magnitude, this is 
not true for «Sg 2 (g), where / + /' is even. This different behavior results from a head-tail 
symmetry which is only slightly broken for our molecules. 

Our main concern has been the investigation of the time-dependent correlators (collective 
and self part) and a test of the predictions of mode coupling theory (MCT). This has been 
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restricted to the diagonal correlators (I = I'). As a by-product we have found that the q- 
and (I, m)-dependence of the self-correlators m (q, t) approximately factorizes, which was 
demonstrated for I = 1, 2 and for q up to 10.6. The reason for this factorization is based on 
a faster relaxation of the ODOF, compared to that of TDOF. 

Concerning the MCT predictions, we first studied the existence of a single transition 
temperature T c . For the (q, I, m)-dependent a- relaxation times Tj mi9 (T) we have found that 
they can be fitted with a power law (1-10) with T c = 0.475 ± 0.01. Thus from the numerous 
correlators we have investigated, one unique temperature T c can be located, at which the 
dynamics of TDOF and ODOF crosses over from an ergodic to a quasi-nonergodic behavior. 
This temperature also agrees with that obtained from the translational diffusion constant 
D(T). This indicates that the TDOF and the ODOF are strongly coupled. Values for 
7 and the corresponding exponent parameter A are given in Table I for the translational 
diffusion constant and a selection of correlators. From this Table we observe that 7 is non- 
universal. Nevertheless there seems to be some systematic behavior. The 7-values for all 
the correlators with / 7^ 1 correspond to A = 0.76 ± 0.03 and are essentially independent 
of q and independent of whether the collective or self correlator is considered. A deviation 
from this value occurs for 7^, the exponent for the diffusion constant, and even a stronger 
one for all correlators with I = 1. A similar discrepancy between 7^ and the exponent for 
the I = relaxation time has been reported before |TjJ, which shows that this prediction of 
MCT seems to be problematic. 

This exceptional role of the (I = l)-correlators is also observed for the first and second 
scaling law of ideal MCT. A consistent picture within ideal MCT emerges for all q,l,m 



with 1^1. The situation is illustrated in Fig. [16] for an exponent parameter A = 0.76. 
There we plot (Su(q,t) — fii m (q))/hfl(q) versus t, which should equal in the first scaling 
regime the critical correlator G(t). All the correlators shown follow the "universal" time- 
dependence of the critical correlator G(t) for A = 0.76. Such a behavior was also found by 
Wahnstrom and Lewis |l^] for a simple model for orthoterphenyl. The time range for which 
the correlators can be fitted by G(t) depends on q, I and m and varies between one and 
a half decade (for C| (t) = S^ m (0,t)) and three decades (for F(q max ,t) = SQ (q max ,t) ). 
Although this time range increases significantly by taking A and the /3-relaxation time scale 
t a as free parameters which seems to yield A — > 1 for I — > 00, we believe that the different 
time ranges relate to the (q, I, m)-dependence of the size of the asymptotic regime. This has 
been demonstrated earlier for the TDOF of supercooled water |Tj| and for the TDOF for 
our molecular system in KKSI. That the asymptotic regime depends on q has recently been 
shown by the analytical calculation of the next order corrections for a system of hard spheres 
|J. We also find that for the correlators with I = I' > (with exception of I = I' = 1) the 
asymptotic regime is largest for q$ ax , the main peak of the static correlator S™(q). This is 
in variance with the result for water [Ol. There it has been found that the corrections are 



smallest for q = qFSDP, where qFSDP is the position of the first sharp diffraction peak and 
not that of the main peak of S(q). This difference probably relates to the different types of 
glass forming liquids. Water is a network former due to covalent bonding mechanism, which 
is absent for our model liquid. The role of this correction to the asymptotic laws is also 
supported by the fact that the (q, I, m)-dependence of the critical nonergodicity parameters, 
shown in Fig. [15, is only consistent with that of fu m {q) obtained from the molecular MCT 
for the present liquid of diatomic molecules, if the next order correction to the von 
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Schweidler law (cf. Eq. (1-9)) is taken into account. 

The result shown in Fig. [16] also demonstrates the validity of the factorization of (q, I, m)- 
and t-dependence of the various correlators on the time scale of t a . For simple liquids, i.e. 
for I = m = 0, this is a prediction of MCT @^J]. There it has been shown that the vertices 
of the mode coupling terms are positive for a simple, one-component liquid, which, however, 
is not true anymore for molecular liquids |§. Since the factorization theorem only requires 
that the largest eigenvalue of a certain stability matrix (see Ref. ||) is non-degenerate, for 
which the positivity of the vertices is sufficient but not necessary, we still believe that this 
non- degeneracy is generic and that therefore the factorization theorem holds for molecular 
liquids as well. In the case that a system exhibits a type-B transition || , this non-degeneracy 
and hence the factorization is guaranteed. 

The exceptional behavior for the correlators with I = 1 has also been observed in the 
susceptibility (cf. Fig. |8|). The position of the minimum between a- and microscopic peak of 
{x u ) i s approximately the same for I = and I = 2, but not for I = 1. For the latter 

it is shifted to higher frequencies by about one order of magnitude. It is interesting that 
this result resembles the experimental results for some glass forming liquids. For instance 
it has been stressed by Cummins et al. [fJTJ , that light scattering data which may include 



contributions from both, I = and I = 2, are consistent with the spectra obtained from 
neutron scattering (which is only 1 = 0), but not with those from dielectric measurements. 



This is nicely demonstrated for glycerol by Lunkenheimer et. al. [22,23]. The situation 



illustrated in Fig. 2 of |23J is exactly what we have found in Fig. |8| for our system. The reader 
should also note that even the relative weight between the intensity of a- and microscopic 
peaks has the same qualitative behavior in both cases, i.e. it is significantly larger for 
I = 1 than for I = and 1 = 2. A similar result has been recently found from a MD- 
simulation of CKN, where the orientational dynamics (self part) of the NO3 ion was studied 



for I = 1 and I = 2 ||24|| . In that paper, and also for the collective dynamics of dipolar hard 
spheres RH], it has been concluded that the different weights of the a- and microscopic 
peaks relate to the different numerical values for the critical nonergodicity parameters. For 
q = is has been argued that < fjf m (due to q = 0, no m-dependence exists) [^4 . 

Since fu' c ^ m {q = 0) is the a-relaxation strength of the corresponding susceptibility and 
(x^ )[T(<? = 0) fulfills a sum rule (on a logarithmic frequency scale), it becomes obvious 
that the ratio between the a-relaxation strength and the area under the microscopic peak 
is larger for I = 1 than for I = 2. Whether this agreement between the susceptibilities of 
glycerol and that for our diatomic molecular liquid is merely accidental or not, is, however, 
not obvious. One has to keep in mind, (i) that dielectric spectroscopy and light scattering 
measures the collective dynamics and not their self part and (ii) glycerol has a permanent 
dipolar moment, in contrast to our diatomic molecules. How far the dipolar interaction 
would change our MD-results is not clear. In addition, we believe that the special role of 
/ = 1 relates to the 180°-jumps of the molecules Whether these jumps exist for glycerol 
also and whether they really cause a shift of the minimum is uncertain. 

To summarize, we may say that the results obtained in Refs. |],|]] and in the present 
paper are consistent with MCT. There is strong evidence for a single transition temperature, 
as it is predicted from molecular MCT || and for the validity of the two scaling laws, with 
exception of the correlators with 1 = 1. Concerning the second scaling regime we have found 
that the 7-exponent is not universal in agreement with earlier work on binary liquids [|TB 
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but at variance with the MD-simulation for water |14| , |T8|j . It will be a challenge to clarify 
the discrepancy for the 7-values. The critical law, which is part of the first scaling regime 
could not be observed, due to a strong interference with the microscopic dynamics. 
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FIG. 6. Time dependence of C\ (t) for I = 2, 3, . . . , 6 (bold lines). Thin lines: /3-correlator 
with A = 0.76. The circles indicate the position of the inflection point. 
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FIG. 10. Time dependence of the collective correlators Sfj l (q,t) for q = 6.5, T = 0.477 and 
I = 1 (solid lines), I = 2 (dashed lines). 
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q = 2.8, I = 2 (filled circles); q = 7.3, / = 1 (open squares); q = 7.3, I = 2 (open circles). T c = 0.475. 
The bold lines represent power laws and the thin lines are a guide to the eye. 
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FIG. 13. S22{q,t) (bold dashed lines) versus t for q = 2.8 , 6.5 and 7.3 for the lowest tem- 
perature and the critical correlator (thin solid lines) with A = 0.76 and a time scale t a = 69. The 
circles indicate the position of the inflection point. 
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ing nonergodicity parameter /^' m (<?) and subsequently divided by the critical amplitude h™(q). 
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TABLES 





D 




c«" 


Ms) 

°6 


Fs {Qrnax) 




F{flmax) 


7 


2.20 


1.66 


2.42 


2.80 


2.56 


2.47 


2.57 


A 


0.67 


<0.5 


0.73 


0.79 


0.76 


0.74 


0.76 



TABLE I. The 7-exponent and the corresponding exponent parameter A (from (1-7) and (1-11)) 
for the translational diffusion constant D and various correlators. 
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